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APPENDIX D 



// sample size 

II xl is a vector of length n, from 0 to 1 with step l/(n-l) 
// initialize a vector of zeros with the same length of xl 
II in ascending order 

// index of samples 



Input: A continuous variable x of dimension w x 7 

Output: 1 . A flag indicating whether the input vector is exponentially distributed 
H = 1 : yes; H = 0 : no 
2. The mean value meanv and minimum value minv of the sample 

Process: 
n = 5\ 

x7 = [0:l/(«-l):l] 
x2 = zeros(l, n) 
B = sorted(jc) 
idx = m* xl 
idx = ToxmdQdx) 

While (idx(i) == 0) 
idxii) = 1 

End While 
idx{ri) = m 
For i = 1 :ai 

x2(0 = B{iclx{i)) 
End For //x2 is the vector of samples 

minv = x2(l); //first element is the minimum value 

meanv = mean(x2); //mean value of samples 
//log-scale x2 
For i = 1 :n 



// make sure indexes are not out of boimd 
// last sample is the maximum value 



Compute x2{i) = l-e 



ji:2(/)-fninv 
minv-meflwv 



End For // if x2 now is uniform distributed, x is exponential distributed 

//later is the KS test, test whether xl and x2 have the "same" distribution 
max__d = 0 
For / = l :n-l 

If (abs(jc2(/) - jc7(/)) > max_d) 

max_d = abs(x2(i) - x7(/)) 

End If 

If (abs(jc2(0 - Jc7(/+7)) > max_d) 

max_d = abs(x2(0 - x7(z+7)) 

End If 

End For 

If (abs(jc2(/7) - x](n)) > max_d) 
max_d = abs(x2(w) — xl{ri)) 
End If 
en = sqrt(A2) 

prob = probks((en + 0.12 + 0.11/en)*max_d) 
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If(prob >0.3) 

H= 1 
Else 

H = 0 
End If 

Return H, nunv, meanv; 
Sorting is done in ascending order. 
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